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ABSTRACT: The Moon has suffered intense impact bombardment end- 
ing at 3.9 Gyr ago, and this bombardment probably affected all of the in- 
ner Solar System. Basin magnetization signatures and lunar crater size- 
distributions indicate that the last episode of bombardment at about 3.85 
Gyr ago was less extensive than previously thought. We explore the con- 
tribution of the primordial Mars-crosser population to early lunar bombard- 
ment. We find that Mars-crosser population initially decays with a 80-Myr 
half-life, with the long tail of survivors clustering on temporarily non-Mars- 
crossing orbits between 1.8 and 2 AU. These survivors decay with half-life 
of about 600 Myr and are progenitors of the extant Hungaria asteroid group 
in the same region. We estimate the primordial Mars-crosser population 
contained about 0.01-0.02 Earth masses. Such initial population is consis- 
tent with no lunar basins forming after 3.8 Gya and the amount of mass in 
the Hungaria group. As they survive longer and in greater numbers than 
other primordial populations, Mars-crossers are the best candidate for form- 
ing the majority of lunar craters and basins, including most of the Nectarian 
system. However, this remnant population cannot produce Imbrium and 
Orientale basins, which formed too late and are too large to be part of a 
smooth bombardment. We propose that the Imbrian basins and craters 
formed in a discrete event, consistent with the basin magnetization signa- 
tures and crater size-distributions. This late "impactor shower" would be 
triggered by a collisional disrup tion of a Vesta-s ized body from this primor- 



dial Mars-crossing population ( iWetherill 



19751 ) that was still comparable 
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to the present-day asteroid belt a 3.9 Gya. This tidal disruption lead to a 
short-lived spike in bombardment by non-chondritic impactors with a non- 
asteroidal size-frequency distribution, in agreement with available evidence. 
This body ("Wetherill's object") also uniquely matches the constraints for 
the parent body of mesosiderite meteorites. We propose that the present-day 
sources of mesosiderites are multi-km-sized asteroids residing in the Hungaria 
group, that have been implanted there soon after the original disruption of 
their parent 3.9 Gyr ago. 

Key words: Moon; Moon, surface; cratering; asteroids, dynamics; mete- 
orites 
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1 Introduction 



Small bodies in the Solar System are usually divided into stable populations 
(like the main asteroid belt, Jupiter Trojans or the Kuiper belt) and the 
unstable, planet-crossing populations (like the near-Earth asteroids, comets 
or the Centaurs). Stable populations are seen as fossilized remnants from 
the formation of the Solar System, while the planet- crossing small bodies 
are short-lived and are constantly replenished from these stable reservoirs. 
There are some exceptions to this classification, like the Trans-Neptunian 
Scattered Disk, whi ch is very long-lived but ey olves appreciably over the age 



of the Solar System (IVolk and Malhotra 



20081 ) . The small Hungaria asteroid 



group, on low-eccentricity, high-inclination orbits j ust beyond Mars, has als o 



2010), 



been suggested as a constantly eroding population (IMcEachern et al. 
but the full implications of the Hungarias' long-term decline have not been 
explored yet. 

The issue of small body stability and instability is important for deter- 
mining the impact history of planets and satellites. The Moon is the best 
studied body in the Solar System after Earth, and reveals a record of the sys- 
tem's early history that has since been erased on our planet. The lunar sam- 
ples returned by the Apollo missions firmly established t hat the Moon wa s 



subject to intense bombardment at 3.85 Gyr ago (Gya) (jTera et al 



1974h 



The nature of lunar bombardment before that time is controversial, as is the 
source of the impactors at 3.85 Gya. The terms "Lunar cataclysm" (LC) and 
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"Late Heavy Bombardment" (LHB) have been introduced to describe this 
bombardment, and the two terms are now used almost interchangeably. A 
non-primordial source of the lunar cataclysm at 3.85 Gya is based on two lines 
of argument. One argument is that the decay of lunar bombardment at the 
time of the large Imbrium impact at about 3.85 Gy a is too rapid to be con - 



Chapman et al. 



sistent with a slowly decaying primordial population 
Independently, dynamical calculations show that the primordia . 



(120071 ). 



ies should have been all but extinct by 3.85 Gya (IBottke et al. 



small bod- 



20071 ). This 



leads to a conclusion that the lunar cataclysm impactor population must 
have become unstable after a period of stability. 

The only plausible mechanism for a delayed instability of a numerous 
small body population would be a late rearrangements of the pl anets. This 



could be a local event, involving an additional terrestrial planet ( IChamberd . 



20071 ) , or global, affecting t 



tern (ILevison et al. 



2001 



re giant planet s and, i ndirectly, the rest of the sys- 



Tsiganis et al 



2005 



Gomes et al 



2005|). While 



there is independent evidence for an unstable 



histo ry (jMorbidelli et al 



2005 



e pisode in the Solar System's 



Nesvorny et al. 



2007 



Morbidelli et al. 



2009, 



20101 ) . there is no independent evidence about the timing of t his instability; 



it is m ost often assumed to coincide with the lunar cataclysm (IGomes et al. 



2005|). 



In 



Cuk et al. 



(120 lOf ) we argue that the record of the cataclysm does not 
match the predictions of the planetary migration-based models. Youngest 
basins (Imbrium and Orientale) seem to have formed much later than the 
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older basins (see section 3), a nd were produced by a non-asteroidal impactor 



population (jCuk et al. 



20101 ). while planetary migration predicts t he late 



2005|). More 



impactors to be derived from main-belt asteroids ( IGomes et al. 
generally, the necessity of a connection between the planetary instability and 
bombardment at 3.85 Gya hinges on the lack of other plausible explanations 
for the late delivery of lunar impactors. In this paper, we will show that 
there should still be a significant number of primordial Mars-crossers present 
at 3.85 Gya, so that this population can indirectly trigger a relatively short- 
lived lunar cataclysm, in full agreement with the available evidence. 



2 Lunar Chronology: Basin Magnetization 



The Moon, unlike Earth, currently does not have a global magnetic field. 
However, lunar surface does exhibit magnetization and these local variations 
in the magnitude and direction of local crustal fields have been mappec 



by instruments flo wn on Lunar Prospector spacecraft (IHalekas et al. 



2001 



Hood et al. 



200ll ). The correlation between the crustal magnetization and 



other surface fe atures (basins in part icular) has been controversial and is cer- 



tainly complex ( IHalekas et al. 



actually be antipodal to basins ([Mitchell et al.l. 



2003) . Some of the ma gnetic anomalies may 



200 8J) and there fore likely 



formed by a transient plasma cloud ( jHood and Artemieval . 120081 . this way 



younger impacts ca n even re-magnet i ze old er basins) , but this is not accepted 



by all researchers (IWieczorek et al 



20111 ) . Here we are interested in lunar 



basin magnetization because it can offer us an alternative way of classifying 



basins that can help determine the chronology of their formation. Magnetic 
anomalies inside basins are generally weaker than those outside basins, but 
are more likely to be of thermoremanent origin and therefore record a past 
ste ady magnetizing fie ld (Lon Hood, personal communication). 



MohitetaL 



( 120081 ) find that the Imbrian basins (as well as the late Nec- 
tarian basin Hertzsprung) have very different magnetization properties from 
older basins . Whi le the sources of ancient lunar magnetism are unknown, 



Mohit et al. 



( 120081 ) pointed out that t here appears to 



De a correlation between 



the basins ' remanent magnetiz ation ([Mitchell et al 



age (cf. 



Halekas et al. 



20081 ) and their relative 



20031 ) . Most interestingly, only the youngest few 
basins have weak or no remanent magnetization. In contrast, some Nectar- 
ian basins, especially Crisium, have magnetic sign atures t hat a re fully consis- 



tent with a global dynamo-generate d dipole field (IHoodl. 



2011 



analyses of Lunar Prospector data (IPurucker and Nicholas 



. More recent 



20101) produced 



better maps that strengthen the dichotomy between the stronger magnetic 
anomalies within most Nectarian basins and much weak er magnetiza t ion o f 



Hertzsprung and Imbrian basins. Still, the conclusio ns of 



are b ased on spherical harmonic model maps (such as 



Mohit ct al 



(2008) 



201lh . 



Purucker and Nicholas . 



2010 ) and will need to be confirmed using direct mapping methods (jHood 



The simplest explanation for these differences in basin magnetization is 
that the lunar magnetic field was absent at the time of formation of basins 
younger than Hertzsprung, but present before that. This hypothesis has 
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profound implications for basin formation chronology. If there was no global 
magnetic field at the time of the formation of Imbrium, all other basins 
that formed during the same short-lived lunar cataclysm spike should record 
the same (weak or zero) remanent magnetization. Therefore, only Imbrium, 
Orientale, Schrodinger, Hertzsprung and (possibly) Mendeleev basins have 
magnetic signatures and relative ages consistent with a nearly co-eval for- 
mation in a spike. We note that this interpretation of magnetic data does 
not require Hertzsprung and Mendeleev to be result of the cataclysm as the 
field may have died long before the impact spike, and they could have formed 
between the end of the internal field and the lunar cataclysm. 

Our hypothesis has an important implication for the age of the Sereni 



tatis basin. Since Serenitatis basin does record magnetic fields ( jMohit et al. 



2008 



Lon Hood, personal communication), our scenario places its formation 



significantly bef ore the 3.9 Gya age o f Ap ollo 17 samp 



interpretation of 



Haskin et al. 



( 11998( 1 and 



Spudis et al 



es. This is in line with 



( 1201 ll ) who argue that 



those samples originated in the Imbrium event, in co ntrast with the conven 



19870 



tional interpretation that they come from Serenitatis (1 Wilhelms et al. 
Therefore, Imbrium provenance of 3.9 Gya old Apollo 17 samples would sup- 
port our hypothesis, while the 3.9 Gya age for Serenitatis would contradict 
our interpretation of lunar chronology. 

One could argue that largest impacts could have affected the heat flow in 



the lunar mantle and effectively shut down the dynamo, making the en c 



namo synchronous with the impact spike (lArkani-Hamed and Olsonl . 



of dy- 



2010h 
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However, lack of magnetic signature inside of the pre-Imbrium Hertzsprung 
basin indicates that the decay of lunar magnetic field preceded the giant 
Imbrium impact, while magnetic anomalies within the Crisium and Sereni- 
tatis basins rule out shutdown by the those large upper Nectarian impacts. 
Hertzsprung and Mendeleev basins were probably the first to form after the 
end of dynamo, but they are relatively small and therefore poor candidates 
for an end- dynamo event. 

Magnetization of lunar surface sampled by Apollo missions rocks also car- 



ries information about the past lunar magnetic field. 



Lawrence et al 



find that there is no firm evidence of lunar dyn amo in (mostly Imbrian-age 



Garrick-Bethell et al. 



(J2009I) 



Apollo samples, despite previous claims, while 
identify traces of a likely lunar dynamo in 4.2 Gyr-old sample. Some re- 
search ers find evidence for a dynamo as late as 3.7 Gya, recorded in mare 



rocks ( iShea et al. 



20101 ) . but it is not clear if those results could be reconciled 
with Lunar Prospector data. In conclusion, we propose that the lunar basins 
coeval with Imbrium basin most likely all lack strong thermal remanent mag- 
netization, and that this property can be used to deduce the chronology of 
basins not sampled by the Apollo landings. 

3 Lunar Chronology: Impactor Size Distri- 
butions 

Crater counting is the basic technique for determining relative ages of plane- 
tary surfaces, with the older surfaces having accumulated more craters. Us- 
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ing absolute crater counts of superposed craters, basins can be grouped into 
early Nectarian (like Nectaris or Moscoviense), late Nectarian (among others, 
Serenitatis, Crisium, Mendeleev and Hertzsprung) and Imbrian (Imbrium, 



Orientale and Schrodinge r). Given that 



gether in crater densities (jWilhelms et al 



ate-Nectarian basins cluster to- 



19871 ). and the limited resolution 



of crater counting, it is imposs ible to securely 



in a chronological order. While 



Wilhelms et al. 



to be the youngest Nectarian basin, 



put all late Nectarian basins 



Fassett et al. 



987) consider Hertzsprung 



(1201 If ) identify Mendeleev 



as being the least cratered. Direct superposition of basin materials can help 
but is not available for all basins (Mendeleev and Hertzsprung do not overlap 
with any other late Nectarian basins). However, it may be possible to com- 
pare the size-frequency distributions of of craters that are thought to have 
for med before and a fter the decay of the lunar magnetic field. 



Cuk et al. 



(201 0J) po i nt ou t that most of the morphologically young "class 



1" craters (jStrom et al 



20051 ) must be of lower Imbrian age. If the "class 1" 
crater size-distribution is representative of the Imbrian impactor population, 
then the Imbrian impactors has a SFD very different from the preceding im 



pactor populations rec orded on the heavily cratered highlands (jStrom et al. 



2005 



Cuk et al. 



20101 ). Separate measurements not dependent on morpho- 
logical classification would help confirm or refute this distinction between 
Imbrian and pre-Imbrian impactors. Ideally, t hese would b e counts of all 



craters superposed on an Imbrian surface unit. 



Stroml (119771 ) report counts 



of craters superposed on the Orientale basin ejecta blanket, and state that it 



11 



is consistent with the SFD of cra ters on lunar highlands (later called "popu- 
lation 1" by Istrom et all I2OO5) . Us ing the Kolmogorov-Smirnov (KS) test 



we find that the 



Strom et al 



( 120051 ) highland and class 1 populati ons are 



i ncom patible at 0.1% level, but neither can be distinguished from the 



Strom 



(119771 ) post-Orientale craters at 10% 



about 200 craters identified by 



evel. Unfortunately, it appears that 



Stroml (119771 ) on the Orientale ejecta blanket 



are not sufficient to clarify the Imbrian impactor distribution, especially as 



Orientale basin "population 1 " crater SFD claimec 



dis puted by some r e search ers (IMalhotra and Strom 



by 


St 


rom 


( 


1977) 


2011; 


Cuk et al. 



is still 



Wilhelms et al. 



201 lh 



(119781 ) classified craters based on their stratigraphic re- 



lationships. They find that Imbrian (and younger) craters had a different 
SFD than tho se formed in th e Nectarian period. Table [T] summarizes the 



comparison of 



Wilhelms et al. 



(119781 ) Nectarian and combined Imbrian, Er- 



atosthe nian and Copernica n ("Imbrian and younger") crater size distribu- 



tions to 



Strom et al. 



(120051 ) highland and class 1 craters (for this comparison 
we used only craters with D > 22.6). Using only these overlapping subsets 
of those crater populations reduces statistical significance, but we can still 
conclude, based on KS tests that Imbrian (and younger) craters are incom- 
patible with the highland distribution (at <0.5%) , and compatible with the 
class 1 craters (>10%). On the other hand, Nectarian craters are compati- 
ble to highland population (>10%) and much less similar to the class 1 set 



(at < 2.5%; significance is low as there are only < 200 class 1 craters with 
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D> 22.6km). 



Finally, 



Wilhelms et al. 



(119871 ) report counts of D>20 km craters super- 
posed on a number of large lunar basins, which are used to determine the 
basins' relative chronology. Unfortunately, almost all basins have less than 
hundred superposed craters of that size, making statistically robust SFD 
comparisons impossible. Since the magnetization data imply that Hertzsprung 
and younger basins may have formed in a separate event from the preced- 
ing basins, we combined the crater counts on Hertzsprung, Imbrium, Ori- 
entale and Schrodinger basins ("HIOS") and plotted them in Fig. [TJalong- 
side complete Imbrian (and younger) and Nectarian crater populations from 



Wilhelms et al. 



(119781 ) . While the HIOS craters are incompatible with the 



Nectarian set (<1%; Table [TJ, they are more compa tible with the 



Wilhelms et al. 



mbrian 



(1978) 



and younger craters (2.5% < p < 5%). Additionally, 
Nectarian and Imbrian and younger craters are incompatible with each other 
at 0.1% ( clS £1X6 Nectarian and Imbrian proper), confirming once more that 
they represent different populations. 

Therefore, both morphological and stratigraphic classifications of craters 
agree that there were two different impactor populations in the Nectarian 



and Imbrian periods. Nectarian perio d was dominated 



impactors (to use the nomenclature of 



Strom et al. 



by "population 1" 



20051 ). which also seems 



to dominate the total population of lunar highland craters, and has SFD 
simil ar to that of main belt asteroids in the relevant size range (jStrom et al. 



20051 ). However, Imbrian impactors ("population 2"), also represented by 
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Strom et al.l (120051 ) "class 1" craters, have a different SF D that is less "to p 



20101). 



heavy" and closer to having the cumulative exponent -2 (jCuk et al.l . 
Comparisons of these populations to the "ground truth" of craters that likely 
formed after the decay of lunar magnetic field ("HIOS" set) favor the popu- 
lation 2 impactors over population 1 ones, despite problems with relatively 
small number statistics. 

In summary, both basin magnetization signatures and cratering data 
agree with a transition between two impactor populations at about the strati- 
graphic level of the Hertzsprung basin, implying that there was something dif- 
ferent about the way the Imbrian basins (together with possibly Hertzsprung) 
formed. In the next section we will explore possible dynamical sources of 
these two impactor populations. 

4 Mars-Crosser Dynamics 

Mars-crossers are small bodies that cross the orbit of Mars, but not those of 
Earth or Jupiter. While Amors are technically Mars-crossers, their perihelia 
(q < 1.2 AU) practically couple them to Earth and their dynamics is similar 
to that of other Near-Earth Asteroids. Also, bodies with a > 2.5 AU or 
in resonances with Jupiter evolve much faster than the population we are 
interested in, and will not be discussed here. 

Mars-crossers are interesting because, while unstable, they have lifetimes 
significantly longer than other planet-crossing populations in the inner Solar 
System. This is simply a consequence of Mars's small size, making it the 
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most ineffective among planets at clearing its orbit (ISoter 



20061 ). If we are 



interested in the lunar bombardment at 3.9 Gya, p r imord ial Mars-crossers 



Wetherill 



19751 ). as a significant 



are clearly a population of interest (cf. 
portion of them may have survived for 600 Myr. The most comprehen- 
sive exploration so far of the contrib u tion o f primordial planetesimals to the 
LC/LHB was done by 



Bottke et al. 



(120071 ). who found that the primordial 



inner solar system planetesimals could not have formed the late Imbrium and 
Orientale impact ba sins, which undoubtedly formed after 3.9 Gya. However, 



Bottke et al 



( 120071 ) considered only bodies on high-inclination orbits with 
perihelia below 1.2 AU, excluding Mars-crossers. Therefore, their results 
need to be reevaluated by including a primordial Mars-crossing population. 
We study the fate of Mars crossers by numerically integratin g an initial 



grid o f 1000 particles using a symplectic integrator SWIFT-rmvs3 (ILevison and Duncan 



19941 ). The test particles have been placed in an semimajor-axis and inclina- 
tion grid, uniformly filling the box with sides of 1.75 — 2 AU and — 20°. All 
particles had e = 0.2, resulting in initial perihelia in the 1.4 < q < 1.6 AU 
range. The eight planets were started on their present orbits and with their 
current masses and sizes, and the timestep was 5 days. The integrations were 
completed on Harvard University's "Odyssey" computing cluster. 

Figure 2 shows the distribution of surviving particles after 300 Myr. At 
that time, 87 of the original thousand were still in the system. Most of the 
particles are widely scattered through the Mars-crossing region, while some 



are Earth and Venus-crossers. At this point we cloned the remaining bodies 
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by replacing each of them by 20 new particles with identical orbital elements, 
save for slightly different semimajor axis. We then continue integrating these 
1740 particles until the time of 2 Gyr, when fewer than 30 particles remain, 
and the results cease being statistically significant. Figure 3 shows that 
population at 1 Gyr, when there are still 105 particles left. The surviving 
particles cluster in the region just beyond Mars with inclinations of 15-30° 
and low eccentricities. Figure 4 compares this population to present-day 
Hungaria asteroids. We find that the core of this late surviving Mars-crosser 
population matches the Hungaria group rather well. 

Figure [5] shows the decay of our simulated population over time. The 
decay can best be described as a transition between two exponential decays. 
At first test- particles decay w i th a h alf-life of about 80 Myr. This is similar 



to find ings of 



Morbidelli et al. 



2001 ) and the long tail found by 



Bottkc et al. 



(120071 ). After about 600 Myr a slower decay predominates, with a half- life 



of about 600 Myr. This is in agreement with t he lifetime of large 



(mostly unaffected by r adiation forces 



similar to the results of 



Milani et al. 



foun d by 



McEachern et al. 



Tungarias 



(BOlOh . 



(120 10h and 



While the overall decay of Mars-crossers and their evolution into proto- 
Hungarias appears orderly, it must be noted that Fig. [5] shows the collective 
behavior of particles across a range of simulations. SWIFT-rmvs3 integrator 
treats the test particles as massless, but they do have a subtle effect on 
planetary orbits, as the behavior of test particles triggers timestep reduction 
during planetary encounters. Using different timesteps at different times will 
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make the simulations diverge slowly even if the planetary initial conditions 
were the same. This feature of SWIFT-rmvs3 is in this case desirable, as it 
allows us to show evolution of particles in "parallel universes" and minimizes 
the effects of atypical planetary evolution in any one simulation. However, 
one must keep in mind that the Solar System had a singular history and 
if we were somehow able to repeat our experiment using the system's true 
history, there were likely deviations from the smooth decay shown in Fig 4. 
Most important departures from the average decay rates are likely to happen 
while the long-term maximum eccentricity of Mars is growing or shrinking 



(ILaskar , 



1994 . 



20081 ). which would destabilize previously non-Mars-crossing 
bodies, or decouple some Mars-crossers from the planet, respectively. 

Two other limitations of this experiment need to be noted: we start with 
the current Solar System, and ignore General Relativity. There could be sub- 
stantial differences between Fig 4 and the real decay rate of Mars-cros sers i 



Mars started out o 



2008 



Bottke et al. 



a less (or more) eccentric orbit than it has now fjCuk 



2010af ). The exclusion of Relativi ty mainly causes Mer- 



cury to be less stab le than it is in the real Solar System (jBatygin and Laughlin 



2008 



Laskar 



20081 ). and we had to throw out a number of simulations due to 
the inner planets becoming unstable. In those cases, we re-scale the results 
for the remaining population at the next benchmark, ignoring all losses to 



planetary instability. 
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5 Collisional Evolution of Mars-Crosser Pop- 
ulation 



Our numerical experiment utilizes massless test particles, so it cannot give 
us a direct measure of the amount of mass involved. To estimate the possible 
mass of this Mars-crossing population, we need to compare our population 
of test particles to lunar impactors at some point of time. For an estimate of 
the mass in these particles, we need to know the size-distribution of bodies 
making up this population. Here we will assume that the primordial Mars- 
crossers had the same size distribution as the asteroids in the inner main belt 
(interior to 2.5 AU). In this range, there are about 150 asteroids larger than 
22 km (which would make a D >300 km basin on the Moon), about 50 bodies 
larger than 60 km (that could make a D > 900 km basin) and a single Vesta 
with D >500 km. Since about 7-8% of all particles eventually hit Earth 
(Table [2]), and about 5% of that number should hit the Moon, we estimate 
that one in ab out 300 of our p a rticle s hit the Moon. Since no basins formed 



after 3.8 Gya (IChapman et al. 



20071 ). there should have been no more than 
about 400 22-km-sized bodies left among Mars-crossers at 700 Myr. As there 
were about 200 Mars-crossers left at this point in our simulation, we estimate 
that each of our second-generation particles correspond to two planetesimals 
of about 22 km in diameter. 

One important reality check of this estimate is to compare the present 
Hungarias to our test particles. Starting with about 100 test particles at 



18 



3.5 Gyr ago and a half-life of 600 Myr, we expect about four 22-km Hun- 
garias to survive to the present epoch, ignoring collisions and radiation forces. 
Hungaria region presently contains th e remnants of one b ody of about that 



size that was disrupted about 0.5 Gva (I Warner et al 



other, somewhat older family present ( IWilliamsl . 



20091 ). with a likely an- 



19921 ) . Therefore we seem to 



produce too many Hungarias by less than a factor of two, a relatively small 
disrepancy given the number of assumption s. While ignoring t he Ya rkovsky 



effect for bodies this size is likely justified ( jMcEachern et al. 



20101 ). we do 



not generally expect asteroids smaller tha n 50 km to survive intact against 



20091 ). so we 



collisions over the age of the Solar System (IMinton fc Malhotral . 
need to address the collisional evolution directly. 

Hungarias experience significant bo mbardment by inner main belt aster- 



oids, especially Flora family members (j Warner et al. 



20091 ). so the fact that 



recent families are seen likely means that two 22 km body should have sur- 



vived the col' 
of 



isional attrition almost down to the present. Using approach 



Wetherilll (119671 ). we calculate that the inner main belt asteroids have a 
probability of about 2.7x 10~ 18 km~ 2 yr _1 of striking Hungarias, and that the 
average collisional velocity is about 10 km/s. Therefore a body of only 0.4 
km across should be able to disrupt a 22-km Hungaria. If there are about 
2 x 10 5 such bodies in the inner main belt, we get about 3-Gyr half-life . Flora 



family has a 4 times higher collisional probability ( IWarner et al 



20091 ). while 



containing about 20% of the inner belt's members, increasi ng the number o 



impacts by about a factor of two for the last Gyr or so (jNesvorny et al 
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20021 ). Therefore our rough estimates are in line with the 22-km Hungaria 
population decaying by a factor of two through collisions between the end of 
the Mars-crosser collisional equilibrium and the present epoch. This is sep- 
arate from the 100-fold dynamical depletion that is independent of asteroid 
size. We conclude that our estimate of the Mars-crosser population based 
on the cessation of basin formation at 3.8 Gya is roughly consistent with the 
observed population of Hungarias. 

Extending this estimate to our initial conditions, we find that the orig- 
inal population of Mars-crossers was 400 times larger than our comparison 
population (asteroids with 2 < a < 2.5), and had a mass 60 times that of 
the present asteroid belt (or 3% of Earth's mass). This mass is rather large 
and probably close to the upper limit of what Mars could have scattered, 
but it is still a small fraction of the mass in terrestrial planets. While much 
more numerous than main-belt asteroids, our Mars-crossers are all interior 
to the swee p of the secular resonan ce, which may have depleted the as- 



teroid belt (IMinton fc Malhotra 



2009 



as we explain the lunar cataclysm 



by other means, we assume that the planetary migration and nu% sweeping 



happened early on) . A 



tern planetesimal disk (IHansen 



ternatively, if there was a t r uncat ion in the inner sys 



2009 



Walsh et al. 



20101 ) . we would expect a 



steep decline of the density of planetesimals as we move from the terrestrial 
planet zone to the present asteroid belt. In any case, the survival of present 
Hungaria asteroids require an initial non-Mars-crossing population interior 
to 2 AU ("proto- Hungarias") similar to that of the inner belt, as current 
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20101 ). 



replenishment mechanisms are not very efficient (IMcEachern et al. 
Unless these low-e Hungarias were planted by some unknown mechanism 
that specifically produces low-e (but high-z) orbits, they had to be captured 
from among the Mars-crossers. As we find the capture efficiency to be on the 
order of a percent (Fig EJ), it seems that we cannot avoid having hundreds 
"Vestas" on Mars-crossing orbits early in the history of the solar system. 

Note that this population initially has about 10 4 bodies that can disrupt 
Vesta, but their half-life is only about 80 Myr, so the probability of Vesta 
being disrupted is only 20%, so the survival of Vesta is not by itself in conflict 
with such a large primordial Mars-crosser population. 

There is at least one serious caveat to this estimate. Our initial con- 
ditions are rather arbitrary and we start the integration with the already 
dynamically "hot" population, which may not be correct interior to the i/g 
sweep. While it is likely that a cold population would eventually get heated 
and end up evolving similarly to our sample, orbital distribution of the ini- 
tial population has profound implications for collisional evolution. Since our 
population monotonically decays, most intense collisional grinding happens 
early on when the Mars-crossing region is most densely populated. In the 
following sections we discuss the number of surviving Vesta-sized bodies at 
later epochs, but most bodies of that size would be disrupted early on if the 
initial population was indeed dynamically "hot" like our initial conditions. 
So while our simulation can be used as a proof that a generic Mars-crossing 
population decays with a half-life of 80 Myr and about 1% of it enters the 
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Hungaria region, the simulation plotted in Figs. 2-5 cannot fully correspond 
to the real history of the Solar System. 

A more natural setup would be a significantly "colder" (and probably 
even lar ger) disk beyond Mars th at Mars scatters and may even migrate 
into (cf. iMinton and Levisonl . I201CH ) . This kind of interaction requires using 
massive planetesimals and is outside of the scope of this work. Therefore our 
calculation for what happens before 4.3-4.4 Gya are not to be taken literally, 
and likely only the evolution of Mars crossers after that time is reasonably 
well described by our simulations. This caveat does not affect our conclusion 
that Hungarias are a natural long tail of Mars crossers, and that the extant 
Hungarias indicate a large but reasonable primordial population in the region 
between Mars and the asteroid belt. 

6 Lunar Bombardment by Mars Crossers 

After we have established an estimate of the Mars-crosser flux, we can com- 
pare it to the lunar record. If the last 300 km basin formed no later than 
3.8 Gyr, when did the last 900+ km basin form? Using Fig. 4, the most 
likely time for this would be about 4.0 Gya. One could argue that the 
formation of Imbrium at 3.85-3.87 Gya was a plau sible fluke, but Imbrium 



was shortly followed by 900-km Orientale impact ( jStoffler fc Ryder 



20001 ). 



While we avoid making a case based on small number statistics, we note 
that this top-heavy Imbrian basin distribution inspired George Wetherill to 
look for alternatives to asteroid impacts as a source for the lunar cataclysm 
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(Wetherill 



19751 ) . When we compare our results with those of 



Bottke et al. 



(120071 ). we fi nd that the formatio n of Imbrium and Orientale (within the time 



limits set by 



Bottke et al. 



20071 ) is a 1.5-si gma (<10%) e vent i n our model, 



Bottke et al 



(120071 ) found using 



significantly more likely in our model than 
NEA-like initial conditions. 

There is, however, a strong argument against Imbrium and Orientale sim- 
ply being relatively low-probability events that happened unexpectedly late. 
Apart from the basins, Imbrian system on the Moon contains a large number 



Wilhelms et al 



of sma ller craters. Using numbers of D >20 km craters from 
( 119781 ). we find that at least 20 such craters formed per 10 6 km 2 in the Imbrian 



period (together with three basins), compared to about 50 per 10 6 km 2 (and 
10-12 basins) in the preceding Nectarian period. Therefore the two largest 
Imbrian impactors were accompanied by (at least) a proportionate number 
of small bodies, and could not have been statistical flukes drawn from the 
same Mars-crossing impactor distribution that formed the Nectarian basins. 

It is thus very unlikely that the primordial Mars-crossers could be the 
chief source of Imbrian impactors. Could they be the main source of Nec- 
tarian bombardment? The size-distribution of Nectarian impactors is cer- 
tainly in agreement with the size-frequency distr ibution of main-belt as- 



teroids ( IWilhelms et al. 



1978 



Strom et al. 



20051 1 . which is also our as- 



sumed size-distribution of Mars crossers, which were in collisional equilib- 
rium. While we do not know the actual SFD of primordial Mars-crossers, 
there is no reason to believe it would be different from that of neighbor- 
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ing main belt asteroids. And after 4.4-4.3 Gya they would have a signifi- 
cantly higher flux t han impactors derived from Earth-crossing planetesimals 



Bottkc et al 



2005 



(120071) or long-term erosion of the asteroid belt (I Gomes et al 



Minton and Malhotra 



2010 



; assuming an early planetary migration). 
Assuming bombardment by Mars-crosser population shown in Fig. [51 
Nectarian period would have lasted from about 4.2 Gya (when the num- 
ber of Mars-crossers has declined by about 20), until the Imbrium impact 
at 3.87-3.85 Gya. In section 2, we argue that basin magnetization signa- 
tures require a substantial gap between the formation of most Nectarian 
basins and Imbrian ones. If the last pre-gap D >900 km basin was Seren- 
itatis at about 4.0 Gya (see above), this would leave about 150 Myr for 
the decay of the lunar magnetic field (which would have already been rel- 
atively weak at the time of Serenitatis impact). This timing of Nectarian 
basin formation is consistent with the existence of lunar dynamo at 4.2 Gya 



(IGarrick-Bethell et al 



namo (jDwyer and Stevenson 



2009) . and the hypothesis of a nuta ti on-driven dy- 



2005 



Meyer and Wisdom 



20111 ) 



Meyer and Wisdom 



( 1201 ll ) find that the strong dynamo was most likely when the Moon was at a 
distance of 34 Earth radii from Earth, which corresponds to about 250 Myr 
after its formation if early Earth's tidal quality factor Q was about 100 (with - 



Bills and Ray 



19991 ) . 



out strong constraints, this value is certainly allowed; 
While detailed chronology (and the issues with dating of Serenitatis) will be 
discussed in Section 8, we will conclude here that the Mars-crossers are the 



likely candidates for Nectarian impactors, as they they plausibly could have 
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had the right size-distribution, they survived later than other primordial im- 
pactor populations and they can satisfy the chronological constrains based 
on basin remanent magnetization. 



7 Mutual Collisions and the Imbrian "Cata- 
clysm" 

Primordial Mars-crosser population was likely insufficient to produce all of 
Imbrian basins and craters. In light of conclusions of Sections 2 and 3, this 
is not surprising, as there we find that there was a change in the impactor 
population between the Nectarian system and the Imbrian system. While 
the Nectarian impactors can plausibly be primordial (and we made a case 
for Mars-crossers), at least some Imbrian impactors had to be released later, 
in an event that we could call a "minimal lunar cataclysm" or "Imbrian 
cataclysm" . 



CuketaL 



(120 Iff ) argue 



What do we know about Imbrian impactors? 
that their size-distribution was not asteroidal and therefore they cannot be a 
product of a gravitational destabilization of an asteroid-like re s ervoir . This 
conclusion is based in part on the findings of 



Wilhelms et al. 



( 119781 ) . who 



found that Imbrian impactors had an excess of small bodi e s com pared to 



Nectarian ones. Regarding their composition, 



Anders et al. 



( 119731 ) thought 



that the Imbrium impa c tor w as closest in trace element si 



(cf. 



Ganapathy et al 



gnatur e 



19721 ) or IIIA iron meteorites. ( 1 James . 



to ty pe IVA 



20021 ) found 



that pre-cataclysm impactors were chondritic, while those at 3.9 Gya were 
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a mix of rarer types, including type IAB irons and pallasites. 



Puchtel et al. 



(120081 ) agrees with the non-chondri tic (possibly iro n meteorite) identification 



for the Imbrium impactor (as does 



Korotev 



19871 ). while finding more chon- 



dritic composition for some other lunar impactors. T 



re rate of impact s after 



Imbrium and Orienta 



Chapman et al. 



e impacts dropped rather fast ( jStoffler fc Ryder 



2000 : 



20071 ). and the rate of this decay is consistent with NEAs 



which are Earth-crossers, Mars crossers or a mixture of those two 



Cuk et al. 



( 120101 ). but not with a very long-lived population (with a half-life > 100 Myr). 

An important factor to consider here are mutual collisions between Mars- 
crossers. We postulate that most of their mass was locked in Vesta-sized 
bodies. A disruption of one of them may be able to both increase the prob- 
ability of large basin formation and modify the size-distribution of smaller 
craters by injecting small impactors with a steep size-distribution into the 
Mars-crosser region. If the Imbrium impactor was about 60-km across, a 
thoroughly-disrupted Vesta sized body could produce hundreds of such im- 
pactors, enough for one of them to strike the Moon. Since there are about 
30 bodies capable of disrupting Vesta at 10 km/s in our reference population 
(i,.e. with D> 90 km), and the population at 550 Myr (one Mars-crosser 
half-life before the Imbrium impact) was about 8 times the reference one 
(with 2% of 400 still surviving), and the mutual collision probability in our 
simulation at 600 Myr was 1.3 x 10 -17 , we get that one of "Vestas" should 
be disrupted every 400 Myr. As our population decays with a half-life of 



about 200 Myr at this point, probability of such disruption at 550 Myr or 
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later is about 25 %. Such an event is not improbable, and naturally explains 
both the new impactor size distribution and increased impact rates during 
the formation of the Imbrian basins. 

How large would the Wetherill's object need to provide the Imbrium and 
Orientale impactors? Since most mass is in the Imbrium impactor and its 
energy was of the order of 10 26 J, and assuming v = 20 km/s, Imbrian im- 
pactor may have had a mass of about 5 x 10 17 kg, and the total mass of 
the Wetherill planetesimal would need to be about 2 x 10 20 kg or larger, 
assuming Mars-crosser-like impact probabilities. The mass of the Orientale 
impactor does not change this estimate much (especially as the largest im- 
pacts suffer from small- number statistics). Smaller impactors can safely be 
ignored, as about 700 Imbrian impactors forming craters with D > 20 km 



likely carried only about 10 16 kg. Imbrian craters in the 12 
have the cumulative size distribution close to D~ 2 ( ICuk et al. 



00 k m range 



20101), imply- 



ing only a few times larger mass in impactors making 100 km craters. Note 
that the extrapolation of this size distribution would predict several 300 km 
basins and no D > 900 basins, implying that basin-forming impactors did 
not follow this power law, and there may have been a "bump" at the size 
of Imbrium impactor (Fig. [6]). If the Imbrium impactor was indeed metal- 
lic, then it seems that even the core of the Wetherill's object was disrupted. 
Therefore, if we assume that most of the mass of the Wetherill's object was 
in the Imbrium impactor-sized fragments, we arrive to a total mass estimate 
close to the mass of Vesta (which is 3 x 10 20 kg). 
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8 Connection to Mesosiderites? 



While the the disruption of the Wetherill's object would explain the pecu- 
liarities of the lunar bombardment in the Imbrian era, additional evidence 
of its late disruption would be invaluable. Apart from the lunar cratering 
record, the only other current window in the absolute chronology of small- 
body populations is offered by meteorites. We have already described the 
constrains on the Wetherill's object from the Imbrian cataclysm: it should be 
about the size of Vesta (D ~ 500 km), formed in an environment similar to 
Vesta's, differentiated, and disrupted at about 3.9 Gya. We propose that the 
mesosiderite group of stony-iron meteorites satisfies all of these constrains. In 
addition to the similarity in inferred properties between mesosiderite parent 
body (MPB) and the Wetherill's object, our hypothesis would also explain 
the puzzling absence of the MPB from the current asteroid belt. 

The mesoside rites are made of about equal mix of iron and silicates 



((Hutchison 



2004]). The silicate component resembles eucrites, and is gen- 



erally thought to ha ve originally forme d close to the surface of a large dif- 



ferentiated asteroid (IScott et al. 



20011 ) . The metallic component is sim- 



ilar to type IIIAB iron meteorites, and has presumably originated in a 
core of a differentiated body. Available evidence points to the iron-silicate 
mix being established early in the history of the sol ar system, presum - 



ably in a dramatic disruption and re-accretion event ( IScott et al 



2001) 



Apart from their unusual composition, cooling histories of mesosiderites 
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are also peculiar. M etallic component appears to have cooled exceptionally 



slowly, 0.02 K/Myr (lHaack et al.l . 119961 ). This has been taken as a sign that 
mesosiderites were buried deep wit hin a large (D < 400 km) asteroid for hun- 



dreds of Myr ( IPetaev et al. 



19921 ). Ar-Ar da ting of the silicate component 



gives ages that cluste r arou nd 3.94±0.1 Gya (IRubin and Mittlefehldt 



Bogard and Garrison 



19931 : 



19981 ) , but with indications that at least part of the Ar- 
Ar ages comes from cooling rather than a discrete event. Still, slow cooling 
inferred from the metallic component (0.02 K/Myr) may be too slow to ex 



plain the the argon closure for the silicates at 3.94 Gya (IBogard and Garrison . 



19981 ). The cooling of silicate portion at 3.94 Gya may have been greatly 
accelerated by a breakup event, explaining the mismatch in cooling rates. 
However, the lack of strong shoc king in meso s iderit es has been used to argue 



against a collisional disruption (lHaack et al. 



1996), and is used as evidence 



that the late argon dates reflect gradual colling only. We hypothesize that 
the lack of shocking is due to the already fragmented rubble-pile structure 
of the body resulting from its ancient disruption in which the differentiated 
layers were mixed. In any case, it is hard to see how mesosiderites could have 
been excavated from deep within the body without some kind of collisional 
disruption. 

Disruption of the MPB on an unstable orbit would also naturally explain 
the absence of large amount of mesosiderite parent body material in the main 
asteroid belt. Even if the mesosiderite Ar-Ar ages are due to cooling and had 
nothing to do with disruption of the MPB, mesosiderites had to be somehow 
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excavated from great depths. Is is hard to see how the MPB would not be 
completely disrupted in the process. Even if it re-accretes, it would have a 
large family which should still be observable, due to its unusual composition, 
in addition to expected orbital clustering. If the mesosiderites are excavated 
from their parent body later than the 3.8 Gya (by which time the asteroid 
belt must have reached its present state), there is no dynaical way to dis- 
perse or hide such a family. This lack of any fa mily is a serious ar gument 



against metal-rich asteroid Psyche being the MPB (IDavis et al. 



19991 ) . Rela- 



tively small size makes the asteroid Maria an unlikely can didate for the MPB 



despi te some spectroscopic similarities to mesosiderites ( IFieber- Beyer et al. 



201ll ). More broadly, if the mesosiderites were exavated from deep within 



the MPB, large amount of olivine mantle material must have been exca- 
vated, too. While the mesosiderite breccias managed to avoid containing 
large amounts of olivine (possibly due to the elements o f the molten core p ref- 



Scott et al 



20011 ). it 



erentially attaching to relatively cool crustal materials; 
is all but impossible to have the whole re-accumulated MPB be olivine-free. 
It is widely agreed that the mes osiderites cannot b e typical of the their par- 



ent body's overall composition (IScott et al. 



200l[ ). so we would expect the 



surviving MPB to be a giant, re-accumulated, olivine-rich asteroid, which is 
clearly not observed. This need to destroy and re-accumulate the MPB at 
least once and possibly twice makes Vesta the unlikely candidate for MPB, 
despite the simila rity in the isotopic signatures of MPB and HED meteorites 



(jHutchisonl . 



20041 ). The isotopic similarity between some IIIAB irons, HED 
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meteorites and mesosiderites more likely implies some amount of uniformity 
in isotopic composition between large differentiated asteroids that were abun- 
dant around 2 AU. In the past, main-group pallasites were also considered 
isotopically indistinguishable from HEDs and mesosiderites, but better tech- 



niques c an now distinguis 



i their isotopic makeup from that of the other two 



groups ([Greenwood et al. 



20061 ). We note that no known iron meteorites 
have very slow cooling signatures like mesosiderites, implying that no pure 
iron fragments from the MPB core survived to the present. This is was ei- 
ther due to chance or, more likely, the core may have been broken-up into 
a smaller number of large objects, none of which survived scattering by the 
planets. 

The fact that the mesosiderite parent body was thoroughly disrupted and 
reaccreted early in its history, likely through a giant low-velocity collision, is 
in line with our concept of a very massive primordial Mars-crosser population, 
if not with out initial conditions (where e = 0.2 for all particles). As we noted 
before, our initial conditions are artificial and it is likely that our Mars- 
crossers-to-be were originally placed on low-eccentricity orbits just beyond 
that of Mars. In the massive population we are imagining, Vesta-sized bodies 
should have been disrupted about once on average, most likely early on. 
However, only if the population was "cold" would Vesta-sized objects be likely 
to reaccrete (producing the mesosiderite mixture), and survive late enough 
to become the cause of the lunar cataclysm. Th e dynamical heating of this 



population by Mars may have involved migration (jMinton and Levison 



2010) 
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and it is clear that our test-particle simulations cannot be taken literally in 
the first few hundred Myr of the Solar System's history. However, they 
prove that after being dynamically heated, Mars-crossers naturally leave a 
Hungaria-like tail. 

9 Survival of Mesosiderites 

A collisional breakup of a Mars-crossing MPB solves the problem of its ab- 
sence from the current asteroid belt. However, a new problem appears, re- 
garding the survival of any of its fragments. If the precursors of mesosiderites 
were planet-crossing, it is surprising that any of them are still around. It 
is clearly im possible that mesos iderites could survive almost 4 Gyr in near- 



Earth space (IBottke et al 



20071 ). However, our own integrations have shown 
that on the order of 1% of Mars-crossers become trapped in the Hungaria 
region, decaying very slowly. Assuming that Hungarias were depleted by 
two orders of magnitude since the fragments' implantation (as about 6.5 
of the 600-Myr Hungaria lifetimes have passed since LHB), we find that 
about 0.01% of the MPB's fragments should have survived among Hun- 
garias, based on dynamical considerations alone. These would include about 
100 km-sized bodies (of one million original ones) or about one 10-km body. 
Collisional evolution of Hungarias complicates this calculation, as km-sized 
cannot survive to present but new ones could have been made from the 
larger bodies that were disrupted before they had chance to escape. While 
the mesosiderites could not have been favored by the gravitational dynam- 
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ics, it is reasonable to assume that the mesosiderites would be selected for 
by collisional grinding over olivine fragments of the mantle. So as long as 
mesosiderites are mechanically stronger, and one of the D>10 km fragments 
that evolved into Hungarias contained a few percent of the mesosiderite ma- 
terial, the survival of mesosiderites is reasonable. 

There is an additional effect that can account for the survival of mesosiderites. 
If the eccentricity of Mars was higher than the long-term average during the 
MPB disruption, and has subsequently decreased, a larger fraction of MPB 
fragments could have been implanted into the Hungaria region. If the ec- 
centricity of Mars was increasing prior to the lunar cataclysm (i.e. went 
through a maximum around 4 Gya), it could explain both the relatively late 
timing of the cataclysm and help the survival of mesosiderites. An increas- 
ing eccentricity of Mars could somewhat change the decay rate of impactors 
showing Fig 4. If the eccentricity of Mars decreases after some of the MPB 
fragments have been injected into the Hungaria region, large fraction of those 
fragments could be basically stable on the age of the Solar System, rather 
than decaying with a 600 Myr half-life expected from the " average" secular 
histories. In any case, MPB fragments are unlikely to make up more than a 
small minority of all Hungaria asteroids. This estimate is based on the fact 
that the cataclysm should have produced about three times as many as many 
bodies as there were left from formation at this epoch, but that only about 
few percent of these MPB fragments could have been expected to decouple 
from Mars and become proto-Hungarias (which were already a large part of 



33 



the background population). In addition to the varying eccentricity of Mars, 



non-gravitational forces could have also contribute d to the survival o : 



McEachern et al. 



MPB 



2010|). 



fragments ( through the " adoption" mechanism of 
However, this mechanism works best for km-sized pieces and does not seem 
efficient enough to dominate over purely gravitational processes. 

We have assumed that the long-term refuge of mesosiderites would be the 
Hungaria region, as we have previously found that it is a natural storage for 



remnants of planetary accretion. However, th ere are other 



Bottke et al. 



possib ilities for 



( 120071 ) indicate 



long-term survival of small sturdy fragments, 
that the high-z orbits in the innermost asteroid belt may serve as long-term 
refuge to a very small fraction of bodies from the inner Solar System. As 
about a million small km-sized bodies would have been injected into Mars- 
crossing orbits by the disruption of the Wetherill's object, some of them 
may have found very low-probability dynamical pathways to long-term sta- 
bility. One remote possibility is that some of the martian Trojans may have 



been captured at this time, as they may 



( iConnors et al. 



2005 



Scholl et al 



3e long-lived but not strictly stable 



20051 ). One advantage the martian co- 



orbitals would have over Hune long-term survivors is a much more 

benign collisional environment on planet-crossing orbits. While we think that 
Hungarias are the most likely current source of mesosider ites, further work on 



20071 ) 



the dynamics and surface composition of Mars Trojans (IRivkin et al. 
should be able to test the possibility that at least some of them originated 
in the Imbrian cataclysm. 
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10 Discussion 



Here we will detail a new scenario of the chronology and sources of early lunar 
bombardment. This scenario is tentative and will need further improvements 
if it is to explain the full complexity of the lunar record. However, we believe 
that by proposing a number of testable predictions we can help advance lunar 
science. 

On the basis of magnetic and cratering data, we conclude that there 
were two populations of ancient lunar impactors, a pre-Imbrian one and an 
Imbrian one. Thes e groups are largely equivalent to Populations I and II from 



Strom et al. 



( 120051 ). although there is some disagreement about the nature of 
Population II. Here we find that the most likely source of Population I were 
the primordial Mars crossers, which would have decayed with the half-life 
of 80 Myr and may have been able in producing D > 900 km basins until 
4.0 Gya. The size-distribution of the Population I impactors was very similar 
to present-day main belt asteroids (MBAs), but there is no reason to think 
that other (now extinct) inner solar system small body populations did not 
have the sam e SFD. It is reasonable to exp ect that MBAs, Mars- crossers, 



the "E-beltM 



Cuk. 



2008 



formation (IBottke et al 



Bott 



te et al. 



2010aJ) and the leftovers from Earth's 



20071 ) would all have had similar size-distributions. 
The above discussion implies that the size-distribution and composition 
of Population I impactors are unlikely to be specific to their source region. 



However, chronology of their impacts should be indicative of their dynamical 
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evolution. If the Nectaris basin is older than about 4.2 Gyr, it could plausibly 
be the product of a primordial impactor population, as proposed here. If, on 



the other hand, it is yo unger than 4 



a cataclysm-type event (IBottke et al 



Gyr ago, it is likely to be a part of 



2010al ). The absolute age of Nectaris 



basin may be impossible to determine from t he existing Apo ll o 16 samples, 



which appear domina t ed by Imbrium ejecta 



Korotev 



(1998) 



WarrerJ (120031 ): 



Norman et al. 



(11983 ); 



Haskin et al. 



(I2010h . The 3.9 Gya Apollo 17 ages 



attributed to Serenitatis may actually reflect the formation of Imbrium , with 



Serenitatis being much older (IHaskin et al. 



1998 



Spudis et al 



20111 ) . Un- 



fortunately, the factors that complicate dating of the nearside Nectarian 
basins may similarly affect future samples from other locations on the Moon. 
Moscoviense may be the only lower Nectarian basin sufficiently distant from 
both Imbrium and Orientale impacts that it was likely not contaminated 
much by their ejecta, and may be (in theory) the best candidate for sam- 
ple return among Nectarian basins. Dating of the largest and oldest South 
Pole-Aitken (SPA) basin would put a upper limit on the age of l unar surface 



featu res and confirm or rule out the total resurfacing scenario (jTera et al. 



19741 ). but would not give us any information about the chronology of much 



l ater impactors. An a ge of about 4.4 Gyr for SPA would allow both "medium" 



(IBottke et al. 



2010al ) and minimal (as proposed here) versions of the lunar 
cataclysm, so at least one Nectarian basin would need to be dated in order 
to distinguish between these two tentative chronologies. 



In contrast, the timing of (lower Imbrian) Population II impactors is 
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well established, as the Imbrium basin formed at 3.85-3.87 Gya, and the 
most Imbrian craters formed shortly afterwards. Since the size distribution 



of Imbrian impactors 



(jWilhelms et al 



1978 



appears to 



3e different from that of Nectarian ones 



Cuk et al 



20101 ) . and different from asteroids, their 



source is unknown. We pr opose a collision al disruption of a Vesta-sized 



Mars-crosser as their source (jWetherilll . 



19751 ) on the basis of t heoretical con 



siderations, tentative composition of the Imbrium impactor ( jPuchtel et al. 



20081 ). and possible connection to mesosiderites. 

A minimal cataclysm from a Mars-crossing collisional disruption also af- 
fected Mars, but the number of superposed craters indicate that the largest 



impact feat ures on Mars 



cataclysm (ITanaka 



1986 



ike Argyre and Hellas are likely to predate the lunar 



Wernerl . 



20081 ). In any case, amount of mass in 200 



km bodies necessary to form these basins is not consistent with the breakup 
of one Vesta-sized body. Our hypothesis predicts that none of the largest 
basins formed at Mars at 3.85 Gya, and, in general, implies much older ages 
for Noachian featu res than the m odels that assume a larger-scale bombard- 



ment at 3.85 Gya ((Werner 



20081 ). Despite the limited extent of the spike at 
3.9 Gya, Mars would have experienced substantial bombardment prior to 3.6 
Gya due to Mars-crossers and this bombardment was the likely the source 



of sh ocking found in the ancient martian meteorite ALH84001 (I Ash et al. 



19961) . 



Shock ages of HED meteorites and H-chondrites a re also sometimes used 



to constrain the early Solar System bombardment (jBogard and Garrison . 
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20031 ; ISwindle et al.l . I2009f ) . Interestingly, both HED meteorites and H- 
chondrites record a much longer bombardment episode than the Moon, mak- 
i ng it hard to ex plain both records with the same impactor population 



(jHartmannl . 



20031 ). We t hink that it is si 



chondrite parent bo dies (IMcSween et al. 



gnificant that the HED and H 



2010 



Bottke et al. 



Gaffey and Gilbert . 



1998 : 



2010U Vesta, and, most likely, Hebe ) are both located in the 
inner main belt and would be likely to suffer significant number of collisions 
with late surviving Mars-crossers/proto-Hungarias. The shock resetting of 
meteorite Ar-Ar ages appears to be dependen t both on the size of t he body, 



parameters of impacts and the family history ( INesvorny et al. 



2009] ) . There- 



fore, while the number of impacts by Mars-crossers onto Vesta and Hebe from 
3.5 to 4.1 Gya did not dwarf their collisions with other asteroids since then 
(H-chondrites also show more recent Ar-Ar ages, while HEDs do not), it is 
likely that the shock ages of HED meteorites and H-chondrites do have a 
connection to Mars-crossers. Even more than meteorite shoc k ages, impact 



2000 



20051) seem a t odds 



events recorded in lunar meteorites (jCohen et al. 
the lunar impact hist ory inferred from Apollo samples (lHartmann 



2003 



Hartmann et al. 



20071 ). Events postdating 3 Gya are recorded in lunar me- 
teorites, long after the heavy bombardment ended. Our hypothesis does not 
offer an explanation for this discrepancy, and we suspect that the solution 
lies in the physics of lunar regolith rather than the dynamics of impactors. 

How can our hypotheses be tested? First, new crater counts (down 
to about D = 8 km) on all non-magnetized basins (Imbrium, Orientale, 
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Schrodinger and Hertzsprung) would confirm or falsify the existence of a 
non-asteroidal Imbrian size distribution, without having to resort to strati- 
graphic or morphological identification of craters. Otherwise, the study of 
mesosiderites and Hungarias could probably test their relationship with the 
lunar cataclysm faster than new lunar samples can be collected. Our hy- 
pothesis predicts that the mesosiderite parent body (or bodies) are likely 
located among Hungarias (although innermost main belt and Mars Trojans 
are also plausible refuges). Remote sensing of asteroid composition has ad- 
vanced considerably in the last few decades, and it is possible that we may 
be able to identify the intermediate mesosiderite parent body, if it is big 
and bright enough (small Hungarias are brighter than small MB As). Radar 
observations may also in principle be able to break the deg eneracy among 



X-type asteroids among Hungarias (jOckert-Bell et al. 



2010h . Additionally, 



if the Dawn spacecraft identifies a mesosiderite unit exposed on Vesta (or if 
mesosiderites are traced to any other large asteroid deep in the main belt), 
our hypothesis about the connection between the mesosiderites and the cat- 
aclysm would be falsified. While the collisional disruption hypo thesis can 



be formulated independently from mesosiderites (jWetherilll . 



19751 ) . the case 



becomes weaker once mesosiderites are excluded from consideration. 



11 Summary 



In this work we propose several hypotheses regarding the lunar cratering 
record and the lunar cataclysm: 
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1. We interpret the magnetization signatures of lunar basins to indicate 
that there was a chronological gap between the formation of most Nectarian 
basins and Imbrium basin. This interpretation hinges on 3.9-Gya Apollo 
sample dates attribut ed to Nectaris and Serenitatis actually reflecting the 



age of Imbrium basin (IHaskin et al. 



1998|). 



2. Size-frequency distributions of lunar craters agree with a two-stage 



bombardment, with most Nectarian craters having an asteroid li 



ce size dis- 



tribution, and the Imbrian ones containing more small impactors (jCuk et al. 



2010). 



3. We find that primordial Mars-crossers (with original a < 2 AU) should 
have a 80 Myr half-life and could have produced more observable lunar craters 
than any other primordial impactor population. 



4. We propose 



Moon (IStrom et al. 



t hat N ectarian and older "Population I" craters on the 



20051 ) were overwhelmingly produced by Mars-crossers. 
This dynamical model predicts South Pole-Aitken basin being about 4.35 
Gyr old, Nectaris about 4.25-4.2 Gyr old and the last giant Nectarian basin 
(possibly Serenitatis or Crisium) being about 4.05-4.0 Gyr old. 

5. A long-lived (but still eroding) remnant of this population is still 
present as the Hungaria asteroid group. This makes the Hungarias primor- 
dial, but not truly stable over the age of the Solar System. 

6. Imbrian basins (and possibly Hertzsprung) formed about 3.9-3.8 Gya, 
in what is a separate cratering spike. We propose that this spike was caused 



by collisional disruption of a Vesta-sized body ("Wetherill's object") on a 
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Mars-crossing orbit. 

7. Wetherill's object was a late-surviving member of the same Mars- 
crossing population that cratered lunar highlands. 

8. Wetherill's object was also the original parent body of mesosiderite 
meteorites. 

9. Intermediate parent bodies of mesosiderites diffused into the Hungaria 
region from Mars-crossing orbits. They have since suffered dynamical and 
collisional attrition, with the latter favoring metal-rich mesosiderites over 
more numerous silicate fragments of the mantle. 

Hopefully, future studies of the Moon, asteroids and meteorites will pro- 
vide data against which these hypotheses can be tested. 
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Reference and 
crater population 


W78 Imbrian 
& younger 


W78 
Nectarian 


W78 Imbrian and younger 
W78 Imbrian craters only 
S05 class 1 craters (D >22 km) 
S05 highland craters (D >22 km) 
W87 craters on HIOS basins 


100% 

> 10% 

> 10 % 
< 0.5% 
2.5-5% 


< 0.1% 

< 0.1 % 

< 2.5 % 
> 10 % 

< 1 % 



Table 1: Results of Kolmogorov-Smirnov two-population tests for differ- 
ent crater po pulations discussed in this pa per. Abbreviations used in the 

" fl2005[ ): W78 - 



table: S05 - Strom et al 



Wilhelms et ail (1l987h . 



Wilhelms et all (11978( 1: W87 
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Planet / 


2 nd generation 


Earth-crosser 


Mars-crosser 


Outcome 


Mars-Crossers 


break-up 


break-up 


Sun 


70.8 % 


81.6 % 


78.7 % 


Mercury 


3.3 


1.8 % 


0.9 % 


Venus 


10.1 % 


5.6 % 


6.5 % 


Earth 


7.6 % 


6.4 % 


5.3 % 


Mars 


2.5 % 


0.5 % 


2.3 % 


Outer/Ejected 


5.8 % 


4.2 % 


6.4 % 



Table 2: Fates of different populations of test particles. Particles were re- 
moved after hitting the Sun, approaching within one radius from a planet, or 
reaching the heliocentric distance of 100 AU. First column plots outcomes for 
second generation mars-crossing particles that were removed from simulation 
featured in Figs 1-3 between 300 Myr and 1000 Myr. The other two columns 
list the end-fates of fragments from typical disruptions on high-inclination 
Earth-grazing and Mars (but not Earth)-crossing orbits. 
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Figure 1: Comparison between the cumulative distribution of different pop- 
ulations of lunar craters with D >20 km. The ( very similar) so l id an d 
long dashed lines plot the crater distributions from IWilhelms et al.l ( 119781 ) : 
Imbrian (solid), Imbrian and younger (long dashes), and Nectarian (short 
dashes). The intermediate cross-dot line plots the combined distribution 
of craters superposed on Hertzsprung, Imbrium, Orientale and Schrodinger 
basins ("HIOS" basins). 
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Figure 2a: Semimajor axis-eccentricity distribution of the remaining first 
generation test particles at 300 Myr. 
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Figure 2b: Semimajor axis-inclination distribtion of the remaining first gen- 
eration test particles at 300 Myr. 




Figure 3a: Semimajor axis-eccentricity distribtion of the remaining second 
generation test particles at 1 Gyr. 
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Figure 3b: Semimajor axis-inclination distribtion of the remaining second 
generation test particles at 1 Gyr. 
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Figure 4: Comparison between the eccentricity-inclination distribution of test 
particles at 1000 Myr from Figure 2 and Hungaria asteroids (filled squares) 
and 35 brightest Hungarias (open squares). All first-generation particles had 
initial conditions the dashed line (0 < % < 20°, e = 0.2). 
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Figure 5: Remaining Fraction of our Mars-crosser population over 2 Gyr 
(solid). The straight dashed lines trace the initial 80-Myr half-life decay of 
Mars-crossers and the later 600-Myr half-life decay of proto-Hungarias. 
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Figure 6: Rough estimate of the size distribution of fragments from the pos- 
tulated Wetherill's object disruption. We assume that much of the mass 
is in Imbrium impactor-sized pieces (about 60 km in diameter), and with 
fragments below about 10 km having the cumulative distribution close to 
N(> D) ~ D~ 2 . Note that this slope is determined from the "population 2" 
crater SFD on the Moon and may not reflect impactor distribution at any 
one point of time. Since the fragments would be undergoing constant col- 
lisional evolution (and "dilution" by the backround Mars-crossers) between 
the parent body disruption and lunar impact, it is possible that the imme- 
diate post-disruption size-disttribution below 10 km may have been steeper. 
A substantial part of the mass could also have been in few largest pieces but 
we have no way of inferring their existence. 




